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We report results of molecular dynamics simulations characterizing the solvation and interaction 
of small apolar particles such as methane and Xenon in imidazolium-based ionic liquids (ILs). The 
simulations are able to reproduce semi-quantitatively the anomalous temperature dependence of 
the solubility of apolar particles in the infinite dilution regime. We observe that the "solvophobic 
solvation" of small apolar particles in ILs is governed by compensating entropic and enthalpic 
contributions, very much like the hydrophobic hydration of small apolar particles in liquid water. In 
addition, our simulations clearly indicate that the solvent mediated interaction of apolar particles 
dissolved in ILs is similarly driven by compensating enthalpic/entropic contributions, making the 
"solvophobic interaction" thermodynamically analogous to the hydrophobic interaction. 



Ionic liquids (ILs) are a new class of solvents for 
use in environmentally benign industrial processes and 
are seen as alternative to toxic volatile organic com- 
pounds [H, |3j S|- The ionic nature of ILs has impor- 
tant consequences for the structure of the liquid on the 
nanoscopic level 0, S]- Spectroscopic evidence is sug- 
gesting the presence and importance of the formation 
of intermolecular cation/anion hydrogen bonds Fa- 
vorable and specific ion/cation interactions seem to in- 
duce the formation of a persistent anion/cation net- 
work, which has been shown to be quite tolerable to 
adding both polar and apolar particles [7|. Quite re- 
cently, systematic measurements of the infinite dilution 
properties for a number of gases, including methane, car- 
bon dioxide, as well as the noble gases have been re- 
ported [1, [3, [13, HI- The experimental data indicate 
that imidazolium-based ILs of type l-alkyl-3-methyl- 
imidazolium bis(trifluoromethylsulfonyl)imide (denoted 
as [Cnmim] [NTf2]) exhibit an anomalous temperature de- 
pendence of the solubility of apolar compounds showing 
a decreasing solubility with increasing temperature [lo| . 
Moreover, it was observed that the anomalous behavior 
is found to be even strengthened with increasing particle 
size. We have recently developed an improved (nonpo- 
larizable) all-atom forcefield for imidazolium based ILs 
of the type [Cnmim] [NTf2] and have shown that a num- 
ber of thermodynamical and dynamical properties of 
the pure IL could be reproduced almost quantitatively 
[1^ . Here we show that our forcefield is also capa- 
ble of semi-quantitatively describing the solvation be- 
havior of small apolar particles. More importantly, we 
clearly demonstrate that the solvation is characterized by 
an enthalpy-/entropy-compensation-effect, qualitatively 
similar to the behavior of the hydrophobic hydration of 
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FIG. 1: Solubility (here given as inverse Henry's constant for 
the case of infinite dilution with kj^ =exp [/3 /^ei^.Gas] / PihRT 
[ij) of Methane and Xenon in [Cemim] [NTf2] at atmospheric 
pressure conditions. The symbols indicate data obtained from 
MD simulations using our IL-forcefield determined from 
the potential distribution theorem [isll . The experimental 
data is according to Maurer et al. [10l |. 



small apolar particles in liquid water l^, JJ, 15, l3|- In 
addition, we determine for the first time the existence 
of a "solvophobic interaction" of apolar particles in ILs 
which is stabilized by entropic and counter-balanced by 
enthalpic contributions, similar to the hydrophobic inter- 
action of small solutes in water. 

We perform constant pressure (NPT) MD simulations 
of imidazolium based ILs of the type [C2mim][NTf2] and 
[Cgmim] [NTf2] at a pressure of 1 bar over the tempera- 
ture range of 273 K up to 907 K |22| of system-sizes of 173 
ion pairs using our IL forcefield |12l |. Simulations of at 
least 10 ns length are employed for each state-point. The 
simulation conditions 23 1 are similar as in Ref. [l3|- All 
simulations are performed with the Gromacs 3.2 simula- 
tion program 2J| . We study the solvation of small apolar 
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FIG. 2: Simulated excess chemical potential /iei, as well as its enthalpic and entropic contributions /iei =hex — Tse.x of (a) 
Methane and (b) Xenon in [C2mim] [NTf2], as well as (c) Methane in water (all at atmospheric pressure conditions). The 
symbols represent the data obtained from the simulations applying the potential distribution theorem p^^. Enthalpic and 
entropic contributions in a,b were derived from fitting the ^ex{T)-da,ta, to a third order polynome (shown as black solid line 
[l^). The experimental data for the solvation of Methane in water shown in (c) are according to Refs. [20l. [2ll|. 
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FIG. 3: a,c) profile of free energy Afiex{r) = w(r) + v{r) for 
the association of two Methane (a) and Xenon (c) particles 
with increasing temperature. Here v{r) is the intermolecular 
potential, whereas 'w{r) denotes the solvent contribution to 
the free energy profile. The shown temperatures are 273 K, 
303 K, 343 K, 383 K , 473K, and 703 K. The arrow indicates 
increasing temperature. b,d) solvent contribution to the pro- 
file of the free energy, as well as the enthalpic and entropic 
contributions obtained for T — 303 K. Methane (b) . Xenon 
(d). 



particles for the for the case of Methane and Xenon rep- 
resented by single Lennard- Jones spheres with param- 
eters (o'Me.Me ~ 0.3730 nm, ejvle-Me/^ ~ 147.5 K, and 
^Xe-Xe = 0'3975nm, k = 214.7 K) [H applying 

Lorentz-Berthelot mixing rules [26[ for all cross-terms. 
The solvation free energy per particle is given by the 
excess chemical potential iiex- We determined iiex for 
the case of infinite dilution a posteriori from the IL MD- 
trajectories applying Widom's potential distribution the- 
orem with Hex = -kT\ii{Vexp{-P^r))) / (V). 
Here is /3 = 1/kT, V the volume of the simulation box. 



and $(r) is the energy of a randomly inserted (gas) test- 
particle at position r. The brackets (...) indicate iso- 
baric isothermal sampling as well as sampling over many 
different positions r. More details about the calcula- 
tion are available in Ref. [2^. The solubility of the 
gas-particles is given as the inverse Henry's constant 
/c^^ =exp [/3ML,Gas] /PilRT [15]) where .^L.Gas is the 
chemical potential of the solute in the liquid phase and 
p[l is the average number density of the ionic liquid. Note 
that the number density of neutral ion-pairs is used here. 

Figure 1 compares the solubility of Methane and Xenon 
in [C6mim][NTf2] as obtained from our MD simulations 
with solubility data recently published by Maurer et al. 
[lo| . We would like to emphasize that without any re- 
finement of the potential parameters, the solubility-data 
is quite satisfactorily reproduced. For low temperatures 
the agreement is almost quantitatively. In addition, the 
particle-size dependence is well described. We would par- 
ticularly like to stress the fact that both simulation and 
experiment indicate an anomalous temperature depen- 
dence, showing a decreasing solubility with increasing 
temperature, very much resembling the solvation of apo- 
lar gases in water 27|. It is the purpose of this letter to 
demonstrate that the analogy to the "hydrophobic effect" 
is much more deep-rooted thermodynamically. 

In Figure 2 we show the excess chemical potentials 
of Methane and Xenon dissolved in [C2mim] [NTf2]. For 
comparison the excess chemical potential of Methane in 
liquid water is given. The increasing chemical potential 
with increasing temperature equivalent to the anoma- 
lous solvation. In analogy to the procedure used in Ref. 
[2^ the data points shown in Figure 2 were fitted to 
a third order polynome [31 and the corresponding en- 
tropic and enthalpic contributions were determined from 
the temperature dependence of the fitted ^exiT) with 
Sex = -{diJLex/dT)T,p and h^x = A^ea; + Tsex- From 
Figure 2a, b it is evident that for the low temperature 
regime the negative heat of solvation h^x is counter com- 
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pensated by a negative solvation entropy Sex- Moreover, 
the heat of solvation exhibits a positive slope, revealing 
a positive solvation heat capacity contribution of about 
22 JK-^mol"^ for Methane and about 26 JK"Vor^ for 
the temperature interval between 300 K and 400 K. Both, 
the entropy/enthalpy compensation effect, as well as the 
positive heat capacity are also qualitative signatures of 
the hydrophobic hydration of small apolar particles in 
water [3, IS]- However, the solvation heat capacity is 
about a factor of five to six smaller, and the correspond- 
ing solvation entropies are about a factor of three to four 
smaller compared to the solvation in water (compare with 
Figure 2c and data in Ref. [1^). In addition, the maxi- 
mum of HexiT), which is observed in water around 410 K 
to 420 K is shifted to about 900 K to 1000 K for the case 
of [C2mim][NTf2]. The latter values are of course hypo- 
thetical, since the real ILs are not chemically stable under 
those extreme conditions. Our simulations indicate the 
anomalous solvation behavior of ILs is simply stretched 
out on a much broader temperature scale compared to 
water. 

In addition to the solvation behavior, we also deter- 
mine the solvent mediated interaction between two (iden- 
tical) gas particles. Therefore we calculate the pro- 
file of free energy for the association process by apply- 
ing the potential distribution theorem [l§| with w{r) = 

-kT\n{V eM-PHrur2m\ri - r2\ - r)) / {V)-2fi,x. 
Here ^{ri, r*2) is the energy of randomly inserting two gas 
particles, fiex is the excess chemical potential of the in- 
dividual gas particles. The solvent mediated interaction 
w(r) is related to the gas- gas pair distribution function 
g{r) according to —kTlngir) = w{r) + v{r), where v{r) 
is the intermolecular pair potential between the two gas 
particles. w(r ) is also sometimes referred to a "cavity 
potential" [29|. Similar to the excess chemical poten- 
tial, w(r) was calculated a posteriori from stored trajec- 
tory data using a Monte Carlo procedure. The calculated 
profiles of free energy, as well as the corresponding en- 
thalpic and entropic contributions for the association of 
two Methane and Xenon particles are shown in Figure 3. 
The minimum of the profile of free energy shown in Fig- 
ure 3a, c represents the state where two particles are in 
close contact. Differing from the the profile of free energy 
for small hydrophobic particles dissolved in water, there 
is no pronounced second minimum existing here. Hence, 
the presence of a clearly defined solvent separated state 
is missing, which is likely to be related to the much larger 
size of solvent molecules involved here compared to wa- 
ter. However, in parallel to the behavior in water, we do 
observe a deepening of the first minimum with increas- 
ing temperature. When comparing the well depth of the 
minima in Figure 3a, c. Xenon exhibits stronger temper- 
ature effect than Methane. To determine enthalpic and 
entropic contributions, we have fitted each r-histogram- 
point of the set of w{r, r)-histograms to a third order 
polynome with respect to T. The corresponding profiles 
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FIG. 4: Schematic diagram of the "solvophobic association" 
process. The contact configuration is stabilized with increas- 
ing temperatures by minimizing the entropy penalty. 



calculated for T = 303 K are shown in Figure 3b, d. The 
contact-state for both. Methane and Xenon is stabilized 
by the entropy part, whereas the enthalpic part mostly 
destabilizing. A simple explanation is based on the as- 
sumption that solvation enthalpy/entropy is mostly due 
to changes of the solvent in the first solvation shell, which 
is represented by the solvent accessible surface (SAS). In 
the contact state the SAS is minimized, leading to neg- 
ative net entropy and positive net enthalpy for the asso- 
ciation of two particles as depicted in Figure 4. Hence 
the association of apolar particles in ILs is driven by the 
tendency to minimize the solvation entropy-penalty (the 
phrase has been borrowed from Haymet et al. [s^]) sim- 
ilar to what has been found for water Conse- 
quently, the larger temperature dependence of the well 
depth of the profile of free energy observed for Xenon in 
Figure 3,a,c is simply a consequence of its larger solvation 
entropy (compare data in Figures 2a, b). 

Our simulations reveal interesting new insights into 
the solvation and interaction of apolar particles in ionic 
liquids. However, many questions remain unanswered. 
First of all, we have not shown what is exactly causing 
the entropy penalty of the "solovophobic solvation" in 
ILs. It might well be that the effect is similar to the hy- 
drophobic hydration in water, where the entropy penalty 
has been largely attributed to the orientational bias put 
on the water molecules while trying to keep the hydrogen 
bond network around an apolar particle mostly intact 
[3]. Since the formation of intermolecular hydrogen- 
bonds is also an important feature in ILs Q, a similar 
mechanism might also apply here. However, it might also 
well be that the tendency to maximize anion-cation con- 
tacts in ILs is introducing an ordering constraint in the 
solvation shell of an apolar particle and thus causing the 
entropy penalty. In addition, for the case of water it has 
been recognized recently that the solvation of small ap- 
olar particles, which is "entropy dominated" , is different 
than for large scale particles, which is "enthalpy domi- 
nated" . Hence there has to be a crossover, which has been 
placed on the < 1 nm-scale as proposed by the theory of 
Lum et al. 1^, 31 1 . The predicted size of the crossover- 
lengthscales were confirmed recently by computer simu- 
lations of Rajamani et al. 1321 and by simulation based 
scaled-particle theory 33|, |3J]. Given the larger size of 
the IL molecules and considering the importance of main- 
taining anion/cation contacts, the crossover-lengthscale 
might be shifted to larger values for the case of ILs. Con- 
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sidering that the anomalous solvation behavior of gases 
increases with solute molecule-size 
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ubility of proteins in most pure ILs 35| might be a con- 



Ill, the poor sol- 



sequence of a pronounced solvation-entropy effect. 
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